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Abstract 

The ground states of artificial molecules made of two vertically coupled quantum rings are studied 
within the spin density functional theory for systems containing up to 13 electrons. Quantum 
tunneling effects on the electronic structure of the coupled rings are analyzed. For small ring 
radius, our results recover those of coupled quantum dots. For intermediate and large ring radius, 
new phases are found showing the formation of new diatomic artificial ring molecules. Our results 
also show that the tunneling induced phase transitions in the coupled rings occur at much smaller 
tunneling energy as compared to those for coupled quantum dot systems. 
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I. INTRODUCTION 



The realization of electronic structures by design is one of the ultimate goals of nano- 
science. By varying the size and geometry of nanocrystallites, it is possible to tune and 
control the confined quantum states and thus the electronic, magnetic and optical properties. 
A very successful example of such systems are semiconductor quantum dots and molecules. 
These artificial atoms and molecules[l| exhibit shell-filling and new molecular many-body 
states which have been confirmed experimentally. These structures can be used as more 
efficient lasers and as hosts for storing quantum information, i.e. qubit, because of their 
atomic-like properties. 

Quantum rings (QRs) are known for the Aharonov-Bohm effect and its persistent 
current Q where the phase rigidity of the electron wave leads to quantum oscillations in the 
current. Such nanometer-sized rings are artificial benzene analogues. Experimentally there 
have been several approaches to fabricate ring structures rangingfrom top-down approaches, 
e.g. nano-lithography (e.g. atomic force microscope patterning)]^, to bottom-up approaches 
as strain induced self-organization J4] and droplet MBE epitaxy [5]. The bottom-up approach 
has been used recently for nano-structured assembly of quantum ring complexes as double 
rings[6] and vertically coupled ringsj?]. Due to the ring geometry, these ring complexes 
open a new route for measurement of quantum interference effects and for novel many-body 
states. 

In the present work we study two vertically coupled rings which exhibit new molecular 
many-body states, hitherto not found in other systems. The tunability of the ring radius and 
the inter-ring distance leads to a very rich variety of many-electron ground states. In previ- 



ous theoretical studies for single rings, Simonin et al. 



showed that a displaced parabolic 



model is suitable for describing the confinement potential for a self-assembled semiconduc- 
tor QR. Many-body spectra of single QRs have been obtained using exact diagonalization 
techniques^ and spin-density-functional theory 10]. A study of the persistent current in 
QRs for up to 6 electrons showed that the current-spin-density-functional theory provides 
a suitable tool for describing the physics related to the Aharonov-Bohm effect llj . Very 



recently, theoretical 



coupled QRs 



12, 



studies were reported on the single-electron spectrum in two vertically 



131 ] and of few electron eigenstates of concentric double QRs 



Two vertically coupled quantum rings (CQRs) form a new type of artificial molecule 



2 



(AM) where the ring radius together with the inter-ring distance are new tunable parameters 
providing new degrees of freedom to modulate and control the electronic structure of the 
artificial ring shaped "molecules" . In the present work, we apply the well established spin- 
density functional theory to study the ground state configurations or phases of few-electron 
CQRs. Full phase diagrams of these new quantum ring AMs are obtained up to 13 electrons. 



II. THEORETICAL MODEL 



The electrons are confined in a plane where the confinement potential is modeled by a 
displaced parabolic function V(r) = \m*u)Q{r — r ) 2 , where r = (x,y) = (r, 9), loq is the 
confinement frequency and r is the radius of the ring. The tunneling between the two 
stacked identical rings in the z direction occurs through a potential V{z) modeled by two 
coupled symmetric GaAs quantum wells with a finite barrier height Vo = 250 meV and a 
well width of W = 120 A. For these parameters we obtain for the energy splitting A = 22.86 
exp[— d(A )/13.455] meV between the bonding and anti-bonding states when the two wells 
are separated by a distance d. 

We use the Kohn-Sham orbitals VVmo-(r) = exp(— im9)<p nma (r)Z(z) to express the density 
and ground state energy. The eigenfunctions (pnmaif) are expanded in the Fock-Darwin basis 
to solve the Kohn-Sham equation. The contribution from higher excited states is neglected 
because the confinement in the z direction is much stronger than that in the plane. Therefore 
the motion in the z direction may be assumed to be decoupled from the in-plane motion. 
The total density in the rings is p(r) = J2 a J2 n ^ m \4 > nmcr( r )\ 2 where a = (t, I) being the 
z component of the electron spin. Following Ref. 16J], we approximate the density in the z 



direction by S functions at the center of the quantum wells. The exchange-correlation energy 
is treated within the local density approximation using the Tanatar-Ceperley [TZfl functional. 
The ground state energy and the corresponding electronic configuration of the system are 
determined by comparing the total energies of all possible configurations. The ground state 
phases are labeled by three quantum numbers (S z , M z , I z ): total spin S z , total angular 
momentum M z and the isospin quantum number I z , which is the difference between the 
number of electrons in the bonding state and in the antibonding state divided by 2. For the 
CQRs, the ground state phases are found for different inter-ring distances (tunneling energy 
A) and ring radii. In the limits of small and large inter-ring distance d, the results for single 
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QRsjlof are recovered and in the limit of small ring radius (ro — > 0), the previous results for 
two coupled quantum dots (CQDs)[l^, Q are obtained. The energy levels e„ jm of a single 
quantum ring are shown in Fig. 1 as a function of the ring radius (in unit ao = \J Ti / m* ujq) . 
The energy levels are labeled by the radial quantum number n = 0,1,2, and the angular 
quantum number m — 0, ±1, ±2, .... As compared to the quantum dot case, in a ring the 
degeneracy of, e.g. the states e ,±2 and ei )0 (eo,±3 and e^i, etc.) is lifted. With increasing 
ring radius, the energy difference e 10 — £o,±2 increases and tends to a constant Huoq at large 
r . The effect of lifting the degeneracy becomes significant for systems with 9 electrons or 
more when the energy levels e^o and eo,±2 become occupied. Notice that the low lying energy 
levels of the m = orbitals exhibit a minimum as a function of the ring radius. This effect 
can be understood qualitatively by looking at the potential height m*uj\r\j'l in the center 
of the ring that is indicated by the thick dot-dashed curve in Fig. 1. On the left side of this 
curve, the system behaves similar to a quantum dot and to the right of it, a quantum ring 
starts to form. The local minimum in the m = curves occurs at this crossover regime and 
will affect the ground state phases for CQRs with 12 and 13 electrons with intermediate ring 
radius. 

We want to stress that the above theory is based on the envelope function approach 
and the local-spin density approximation (LSDA). These approximations are valid in the 
present study for the coupled quantum rings. For a confinement strength of huj = 5 meV, 
the corresponding typical length scale is ao — 15 nm in the case of GaAs[3|. We show 
results for a diameter up to 2 x 5xao = 150 nm. These scales are still much larger than 
the interatomic distances, therefore the effective mass approximation is valid. Moreover, we 
have chosen for a LSDA-DFT approach to study the many-electron states of the CQRs. This 
approach was used before successfully to describe the energy levels in a single quantum dot 
of comparable size 18]and worked surprisingly well even for systems of 2 and 3 electrons. 
Ferconi and Vignale 18( compared the exact result for the ground state energy of a two 
electron quantum dot with a DFT calculation. The accuracy is found to be better than 
3%. Therefore we believe that claims we make that are based on the LSDA-DFT approach 
are justified. We did not implement a higher functional, because it would not reveal new 
physics. Our confinement potential is anyhow a model system for two coupled rings. 
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III. NUMERICAL RESULTS AND DISCUSSIONS 



The stability region of the different ground states are summarized in phase diagrams in the 
A versus r plane for systems with fixed electron number N (see Figs. 2(a-d)). The insets in 
the figures give the single-particle picture representation for each phase. In the calculations, 
we took the confinement energy huo = 5 meV. Several of the transitions between the phases 
can be understood from a single-particle (SP) picture when the exchange-correlation effects 
are added. The confinement in the z direction splits the single QR levels (Fig. 1) into a set 
of bonding and antibonding levels by A. 

^From such a SP picture, one does not expect configurations different from the CQDs as 
long as the level e 1)0 is not occupied. Indeed, for CQRs with iV < 8 we found qualitatively 
similar ground states as for CQDs. For N = 8, four different phases are found as shown in 
Fig. 2(a). For fixed r , there is a single ring atomic type phase (2, 0, 0) at small A where the 
inter-ring distance is large and the rings are decoupled. For large A the two rings are strongly 
coupled and act as a single one resulting in another atomic type phase (1,0,4). There are 
two molecular type phases: (0,0,2) and (1,2,3) for intermediate inter-ring distance similar 
to those found earlier for CQDs. It is clear that, with decreasing A, more electrons from the 
antibonding state transfer to the bonding state leading to phase transitions. Because the 
SP levels eo, m (for m = 0, ±1, and ±2) decrease with increasing ring radius r and converge 
at large ring radius, the atomic type phase (1,0,4) of strongly coupled QRs is enhanced at 
large tq. To this end, we can also understand that the tunneling effect is enhanced at large 
r . This becomes obvious when we investigate the phase transitions with increasing ro for 
a fixed A. At A = 2 meV (d = 32.8 A), e.g., the transitions (0,0,2) -> (1,2,3) -> (1,0,4) 
corresponds to the isospin transitions I z = 2 — > 3 — > 4 where the maximum value is reached 
for large r . Notice that the magnetic state, i.e. the M z value, can be manipulated by 
varying A and/or r . 

Novel phases appear when N > 9. Let's first study the case of 9 electrons (see Fig. 
2(b). The phase (3/2,0,9/2), for r < 0.89 ao, is consistent with the result obtained for 
CQDs[l6j. A ring atomic type phase (1/2,2,9/2) appears for r > 0.89 a in the strong 
tunneling regime. The emergence of the new phase (1/2,2,9/2) originates from the lifting 
of the degeneracy of the SP levels ei and e ,±2 in the bonding states. When the ring radius 
increases, e^o — eo,±2 increases and when this energy difference overcomes the exchange 
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energy due to Hund's rule, the state eo,±2 becomes energetically more favorable than e^o- 
Consequently, a ground state transition undergoes from (3/2,0,9/2) to (1/2,2,9/2). 

Another new result for the CQRs of 9 electrons is the disappearance of the ground state 
(1/2,2,5/2) for r > 3.04 ao which can be understood as follows: (i) the energy difference 
e o,±i — e o,o i n the antibonding states is small at large ro (it is much smaller than eo,±2 — e o,±i); 
and (ii) the existence of the (3/2, 0, 7/2) phase is largely due to the exchange gain of the two 
aligned spin states in eo,-2 and eo,+2 in the bonding states. For r > 3.04 ao, this exchange 
gain becomes larger than e ,±i — e . With decreasing A, as soon as one of the electrons 
in the eo,±2 levels in the bonding states jumps to the antibonding state, the other electron 
follows leading to the transition from (3/2,0,7/2) directly to the (1/2,1,3/2) phase. This 
is a clear manifestation of the more pronounced effect of the electron-electron interaction in 
CQRs. 

For 11 (and 12) electrons the number of possible ground states jumps to 8 (see Fig. 
2(c)). Among them, the ring atomic type phase (1/2,3, 11/2) and the ring molecular type 
phase (1/2,0,9/2) appear only at finite ring radius. The transition from (1/2,0,11/2) to 
(1/2,3,11/2) at r = 2.05 ao is associated with the transfer of a bonding state electron 
from e 10 to eo,±3- This is almost completely determined by the SP energy crossover of the 
two levels at r = 2.07 ao. We also notice that the molecular type phase (1/2,2,7/2) is 
completely suppressed at large ring radius ( r > 3.17 a ). A similar phenomenon was 
already found for the N = 9 case and is a consequence of the fact that the energy difference 
eo,±i — eo,o in the antibonding state at large r becomes smaller than the exchange gain in 
the (3/2,1,5/2) phase. We also observe a stabilization of the ring molecular type phase 
(1/2,0,9/2) at intermediate ring radius where the inter-ring tunneling energy A leads to a 
crossing of the lowest level eo,o in the antibonding state with the level eo,±3 in the bonding 
state. The lowest antibonding level e ,o is almost flat with increasing r while the level 
eo,±3 decreases rapidly. As a consequence, the phase (1/2,0,9/2) with one electron in the 
antibonding level eo,o is enhanced. 

The number of possible phases increases to 13 for a CQR with 13 electrons (see Fig. 2(d)). 
From left to right on the top of the figure, we see 3 phases: (1/2,3,13/2), (3/2,0,13/2), 
and (1/2, 3, 13/2)* where the star denotes that, in this new ring type phase, the levels e ,±3 
of higher angular momentum are occupied instead of the level ei in the quantum dot type 
phase (1/2,3,13/2). The phase (1/2,3,13/2)* is due to the SP level e ,± 3 which becomes 
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lower in energy than the level e^o in a quantum ring. The intermediate phase (3/2, 0, 13/2) 
is due to many-body exchange-correlation effects. With decreasing A, one of the electrons in 
the bonding state in the above configurations jumps to the lowest antibonding level leading 
to three molecular type phases: (1/2, 0, 11/2), (3/2, 3, 11/2), and (3/2, 0, 11/2). 

With further reducing the inter-ring tunneling energy, transitions occur from (5/2, 0, 5/2) 
to (3/2,2,5/2) at r = 0.99 ao and to (1/2, 1,7/2) with increasing ring radius. These ring 
molecular type phases are a consequence of the lifting of the degeneracy of the levels eo,±2 and 
ei 5 o at finite ring radius and results also in a more stable (3/2, 1, 3/2) phase. The molecular 
type phase (1/2,3,9/2) turns to be unstable for r > 3.89 ao due to a strong exchange- 
correlation gain in its neighbor phase (3/2,0,11/2). For CQRs of 13 electrons, we found 
that the following transitions are invoked by state changes of two electrons: (1/2, 0, 9/2) to 
(5/2,0,5/2), (5/2,0,5/2) to (1/2,1,7/2), and (1/2,1,7/2) to (3/2,0,11/2) where the total 
spin changes with \AS Z \ > 1. 

Experimentally, it is more convenient to tune the number of electrons 191] occupying the 
CQR. Therefore, we calculate the addition energy as a function of the electron number in 
strong, intermediate and weak coupling regimes as shown in Fig. 3. For r$ = 4.0 ao and 
A = 1.6 meV (d = 35.8 A ) in Fig. 3(a), the two rings are strongly coupled acting as a 
single one. For a single ring, the closed shells are formed 10J at N =2, 6, 10, ... resulting 
in peaks in the addition energy. From the inset we notice that I z increases monotonously 
with increasing N because all the electrons are in the bonding states. When we reduce the 
ring radius to r$ = 2.0 ao while keeping the same tunneling energy A = 1.6 meV as shown 
in Fig. 3(b), we find the CQRs entering the molecular type phase for N > 7. The strong 
addition energy peak at iV = 6 indicates that the system transfers from a closed shell ring 
atomic type phase (0,0,3) into a molecular type phase (1/2,0,5/2) for iV = 7 where there 
is one electron in the lowest antibonding level eo,o- From N = 7 to 8, I z changes from 5/2 to 
4 as shown in the inset of Fig. 3(b) with a total isospin change \AI Z \ = 3/2 > 1/2 resulting 
in a strong reduced current peak Q] because of isospin blockade. In this intermediate 
coupling regime, the ground state phase depends strongly on both the ring radius and the 
inter-ring distance as shown in Fig. 2. This is clearly reflected in the addition energy. The 
peak in the addition energy at N = 12 in Fig. 3(b) corresponds to the molecular type phase 
(0, 0, 4) where there are ten paired electrons in the bonding levels and two paired electrons 
in the antibonding level eo,o- 
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When the ring radius is kept as a constant, the CQRs become weakly coupled for large 
inter-ring distance. Such a situation is presented in Fig. 3(c) for r = 2.0 a and d =67.4 
A(A = 0.15 meV) where the CQRs are in the weak-coupling regime. Because the system 
evolutes to two separated rings at d — > oo (A — > 0), one should find closed shells at N =4, 
12, 20, ... for single rings. We see these features in the inset in Fig. 3(c). However, we also 
find a pronounced addition energy peak at iV=8 indicating coupling of two half-filled rings 
with total spin S z =2. The peak at iV=12 corresponds to two closed shell rings of 6 electrons 
in each ring. 

IV. CONCLUSIONS 

In summary, new quantum ring molecular type phases as well as ring atomic type phases 
are predicted for the CQRs which can be probed experimentally through current measure- 
ments. Our results show transitions of these artificial molecular phases from the CQDs 
(r = 0) to the CQRs and also from two strongly coupled quantum ring AMs to their dis- 
sociation in the weak tunneling regime. In comparison to the system of CQDs, quantum 
tunneling effects are enhanced in the present system where the AMs require less tunneling 
energy to entering the molecular type phases. 
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Figures 



FIG. 1: (color on line) Electron energy levels e„ irrt in a single quantum ring as a function of the 
ring radius. The potential height in the center of the ring is indicated by the thick dot-dashed line. 

FIG. 2: (Color on line) The ground-state phase diagrams of the AMs composed by CQRs in the 
plane of the inter-ring tunneling energy (inter-ring distance) vs. the ring radius for (a) N=8, (b) 
N=9, (c) N=ll, and (d) N=13 electrons. Tiojq = 5 meV. The grey (green) region indicates the 
phases which can be found at tq —>0 (ro — > oo). The yellow region shows the phases which exist 
only for intermediate ring radius. The insets indicate the single particle picture representation 
for each phase where the bonding (antibonding) states are represented on the left (right) of the 
vertical dotted line. 

FIG. 3: The addition energy as a function of the number of electrons in the CQRs with (a) ro = 
4.0 a and A = 1.6 meV; (b) r = 2.0 a and A = 1.6 meV; and (c) r = 2.0 a and A = 0.15 
meV. The insets show the corresponding total spin S z , angular momentum M z and isospin I z . 
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